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Abstract 

We take a critical view of the treatment of threshold effects in SUSY spec- 
. trum computations from high-scale input. We discuss the two principal methods of 

I (a) renormalization at a common SUSY scale versus (b) integrating out sparticles 

at their own mass scales. We point out problems in the implementations in public 
spectrum codes, together with suggestions for improvements. In concrete examples, 
we compare results of Isajet7.72 and SPheno2.2.3, and present the improve- 
' ments done in Isajet7.73. We also comment on theoretical uncertainties. Last 

r~| ' but not least, we outline how a consistent multiscale approach may be achieved. 



p i' 

(p : 1 Introduction 



It is often argued [1-4] that, from measurements of supersymmetry (SUSY) at the LHC 
^ ■ and ILC, the model parameters can be extracted precisely enough to test the high-scale 
structure of the theory, in other words to test the boundary conditions of the soft SUSY 
breaking (SSB) terms. This requires relating, with very high precision, sparticle properties 
measured at the TeV energy scale with the Lagrange parameters at the high-energy scale, 
often the GUT scale. 

Of course, the question arises as to what theoretical uncertainties are involved in this 
exercise. Such uncertainties originate from truncating the perturbation series of (a) the 
running of the DR parameters between the electroweak (EW) and the high-energy scale, 
and (b) the relation at the EW scale between DR and on-shell SUSY parameters. They 
have been investigated in [5] by comparing different state-of-the-art spectrum computa- 
tions. Differences at the level of a few per cent have been found, part of which have 
been traced to higher-order loop effects. Since then important improvements have been 
made in all codes. When comparing today the latest versions of Isajet [6], SoftSusy [7], 
SPheno [8] and Suspect [9], the typical spread in the results is < 1%, which is compat- 
ible with the expected precisions at the LHC and getting close to those expected at the 
planned International Linear Collider (ILC) (see [10] for an online comparison). 



1 



One exception is the mass of the hghtest neutrahno. On the one hand, m^o is expected 
to be measured with per-mille precision at the ILC. On the other hand, in mSUGRA with 
a bino-hke LSP, one finds differences of few per cent in the m^o obtained from Isaj et 7 . 72 
as compared to m^o obtained from Sof tSusy 2 . 0, SPheno 2 . 2 . 3 or Suspect 2.3.4 (these 
latter three programs typicaUy agree to ~ 0.5% on m^o). In this context it should also 
be noted that a ~ 1 GeV error in the mass of the LSP can translate into a ~ 10% error 
in the prediction of its relic density [11]. 

The public spectrum codes under consideration all have 2-loop renormalization group 
(RG) running implemented. A major difference between Isajet on the one side and 
Sof tSusy, SPheno and Suspect on the other lies in the treatment of threshold effects for 
computing the sparticle pole masses. This is the topic of this paper. It was the above 
mentioned discrepancy in the prediction of the LSP mass that motivated this study. Our 
analysis, however, applies to the computation of the pole masses of mixing sparticles in 
general. 

The paper is organized as follows. In Sect. 2 we briefly explain the two methods 
used to compute the sparticle pole masses. In Sect. 3 we review these methods for the 
example of the neutralino sector. We compare results of SPheno 2.2.3 and Isajet 7.72, 
and discuss some shortcomings in the computations as well as ways of improving. In 
Sect. 4 we then analyse the scalar top sector. An improved scheme for the computation 
of sparticle masses in Isajet is presented in Sect. 5. The technical difficulties involved 
with a consistent multiscale approach are discussed in more detail in Sect. 6, and finally 
Sect. 7 contains our conclusions. 

2 Two methods 

The RG equations are employed in the DR scheme. The SUSY mass parameters obtained 
from the RG evolution are hence DR running ones. In order to obtain the physical 
sparticle masses one therefore has to perform the proper shifts to the on-shell scheme. 
The public SUSY spectrum codes can in fact be classified by their method of determining 
the sparticle pole masses. 

One approach, adopted by SoftSusy, SPheno and Suspect^, is to assume that the 
full set of MSSM RGEs (gauge and Yukawa couplings plus soft terms) are vahd between 
the scales Mz and Mqut- The gauge and Yukawa couphng boundary conditions are 
stipulated at Q = Mz, while soft term boundary conditions are stipulated at Q = Mgut- 
The RGEs are run itcrativcly between Mz and M^ut until a convergent solution is found. 
The DR parameters arc extracted at a common renormalization scale Q, usually taken to 
heQ — MsusY = ^/^WJ^H^ or ^mi^rrii^. The one-loop (logarithmic and finite) self-energy 
corrections [12] are then added at that scale ^. We will refer to this method as 'common 
scale approach'. It is relatively straightforward and certainly self-consistent, but misses 
two-loop logarithmic contributions between the renormalization scale Q and the actual 
mass scale of the sparticle. Such logarithmic corrections could be relevant in cases where 
the mass scales involved are severely split, such as in focus-point super symmetry. 

The other method, proposed in [14, 15], is to adopt a multiscale effective theory ap- 

^Early versions of Suspect used the step-beta function approach described below; the method of a 

common renormalization scale is used from version 2.0 onwards. 

^For the neutral Higgs masses and the ji parameter, also the two-loop corrections of [13] are included. 



2 



proach, and to functionally integrate out all heavy degrees of freedom at each particle 
threshold. As this is realized by the use of a theta function at each threshold, we call it 
the 'step-beta function approach' (implying continuous matching conditions for the re- 
maining parameters). The program Isajet does adopt a hybrid technique along these 
lines. In Isajet 7.72, the full two-loop MSSM RGEs are employed between Mz and 
Mqutj with the exception that one- loop step-beta functions are adopted for gauge and 
Yukawa couplings. (This approach means that log corrections to gauge and Yukawa cou- 
plings at the scale Q = Mz are not needed, since they are handled by the RG evolution). 
In addition, each SSB parameter rrii is extracted from the RGEs at an energy scale equal 
to rrii = mi{mi), with the exception of parameters involved in the Higgs potential, which 
are all extracted at the common scale Q = MgusY = 

The RG evolution of mass parameters from high to low scales is equivalent to comput- 
ing the logarithmic radiative corrections. The step-beta function approach hence gives a 
leading-log approximation of m, = mi{mi). For sparticles that do not mix, i.e. gluinos 
as well as squarks and sleptons of the first and second generations, this approach di- 
rectly gives the pole masses up to constant terms. For the required level of precision 
these constant terms are important, so they have to be added at the end of the running. 
The situation is more complicated for mixing sparticles, since multiple mass scales can 
be involved in the mixing matrices, and care has to be taken to compute the on-shell 
mass matrices in a consistent way. This will be discussed in more detail in the following 
sections, first for the example of the neutralino sector, and then for the case of scalar 
tops. Moreover, as we will discuss in Sect. El at the two-loop level, each time a sparticle is 
integrated out this implies non-trivial matching conditions for the parameters remaining 
in the effective theory. To our knowledge these matching conditions are not yet known 
and, thus, are not yet taken into account properly. There are also other complications, 
which we will discuss in Sect. El 

3 Neutralino sector 
3.1 Neutralino mass matrix 

At lowest order, the neutralino mass matrix in the basis ipj = {—i^'i ~'^^^yi'Hiyi'H2) 

/Ml -mzSwC/3 mzswsp \ 

M2 mzCwCfi -mzcwsp , . 

—mzSwCp mzCy/Cp — 

\ iTLzSw-Sp -mzCwSfs -/i / 

with sw = sm9w, cw = cos 9w, S/3 = sin/3, C/j = cos/5 and tan/3 = V2/V1. It is 
diagonalized by a unitary mixing matrix N: 

NAiNN"^ = diag(£im^o, e2m^o^ ^s^^xi' ^■^'^x") ' (^) 

where m^o {i = 1,...,4) are the neutralino masses, m^o < .... < m^^o, and Si are their 
signs. The mass eigenstates are = Nijip'j. 

The SUSY parameters taken out of the RG running are DR parameters at a certain 
scale Q. Straightforward diagonalization of Eq. (P) therefore gives the DR running mass 
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Figure 1: Dependence of the masses of Xi (left) and X34 (right) on the renormahzation 
scale Q. The dashed lines are the DR running masses, the full lines are the one-loop 
corrected pole masses. Computed with SPheno 2.2.3. 

eigenvalues. In order to obtain the neutralino pole masses, one has to add self-energy 
corrections 

X— ^ ^ AMNiQ) , (3) 

leading to corrections in the masses, m^o m^o + Am^o, and in the mixing matrix 

A-Z A-i ' 

N —>■ N + AN. Here notice that the on-shell condition = has to be fulfilled for each 
m-j^o separately; so one has to diagonalize Eq. (jHI four times^. The corrections at the one- 
loop level are given in [12]; see also the discussion in [16]. They typically amount to a few 
per cent. Since the shift from the DR to the on-shell sheme, Eq. (jS)), cannot be performed 
to all orders, there will always be a small residual scale dependence of the pole masses. 
This scale dependence is often regarded as an estimate of higher-order corrections. 

In the step-beta function approach, the mass parameters that enter Eq. are Mi (Mi), 
^2(^2), and /i(MgusY)- According to [15], this corresponds to the on-shell at up to 
finite corrections. The on-shell mass matrix at the full one-loop level is then given by 

_y^onshell ^ + AM^"' ■ (4) 

The complications involved in this procedure will be discussed in Sect. 



3.2 Threshold corrections at a common scale 

We first discuss results for the approach of running all SUSY-breaking parameters to a 
common renormahzation scale Q and adding one-loop threshold corrections at this scale. 
As already mentioned, this is the method employed in SoftSusy, SPheno and Suspect. 
We use the mSUGRA benchmark point SPSla with 

mo = 100, mi/2 = 250, = -100, tan/? = 10, /i > 0, = 175 (5) 

■^This actually leads to a (numerically very small) ambiguity in the neutralino mixing matrix. In 
Isaj et 7 . 69-7 . 72, tree-level mixing elements are adopted for cross section calculations, along with loop- 
corrected neutralino masses. In SPheno 2. 2. 3, the mixing matrix obtained for = m?o is used. 
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Q 


Ml 


M2 


A* 


Xi 


X2 


Xa 


X4 


pole 
Xi 


pole 

X9. 


pole 

X?. 


pole 

m^o 

X4 


100 


98.39 


188.9 


352.7 


95.27 


173.3 


359.3 


377.9 


96.82 


179.9 


365.2 


382.2 


200 


99.78 


190.1 


355.0 


96.66 


174.8 


361.4 


379.9 


96.97 


180.3 


365.0 


382.2 


^SUSY 


101.6 


191.7 


357.6 


98.47 


176.6 


363.9 


382.2 


97.11 


180.7 


364.9 


382.2 


1000 


103.1 


193.0 


359.2 


99.95 


178.0 


365.4 


383.5 


97.16 


181.0 


364.6 


381.9 


Table 1: 


Parameters and masses (in 


GeV) for the mSUGRA benchmark point SPSla 


obtained with SPheno 2 


2 . 3 for different renormalizat 


ion scales Q; Msusy = 


484.5 GeV. 


Q 


Ml 


M2 


/i 


X? 


x" 


x!^ 


X4 


pole 

X? 


pole 


pole 

x§ 


pole 

m-o 
xS 


100 


98.55 


188.9 


352.6 


95.25 


173.3 


358.5 


377.2 


97.00 


179.7 


363.4 


382.4 


200 


99.96 


190.2 


354.8 


96.66 


174.7 


360.5 


379.1 


97.18 


180.2 


363.0 


381.9 


^SUSY 


101.5 


191.6 


356.7 


98.38 


176.4 


363.0 


381.3 


97.31 


180.8 


362.8 


381.7 


1000 


103.3 


193.1 


358.9 


99.96 


177.9 


364.4 


382.7 


97.37 


181.2 


362.1 


381.0 



Table 2: Same as Table [H but computed with Suspect 2.3.4; Msusy = 465.4 GeV. 



as an illustrative numerical example. Moreover, we take as{Mz) = 0.1172, and mh{mh) = 
4.214 GeV as SM input values in the MS scheme so as to be able to compare with 
Isajet 7.72. 

The dependence of the DR running (tree-level) neutralino masses and the one-loop 
corrected pole masses on the scale Q is shown in Fig. for Q = 100 GeV to 1 TeV. In 
addition. Table ^ lists the DR mass parameters together with the tree- level and one- loop 
corrected neutralino masses at SPSla for some particular choices of Q. The numbers have 
been obtained with SPheno 2.2.3, which has the complete one-loop corrections of [12] for 
all sparticle masses. As can be seen, going from the DR to the on-shell scheme is a quite 
important correction. The scale dependence of the DR (or tree-level) neutralino masses, 
m?o^, is about 5 GeV between Q = 100 GeV and Q = 1 TeV. That is about 5% for the 

Xi 

LSP mass and about 2-3% for the masses of the heavier neutralinos. For the one-loop 
corrected pole masses, m^o^*^, the scale dependence goes down to the level of a few per 

Xi 

mille. Here notice also that the scale dependence is largest for the wino-like X2- 

The results of Sof tSusy 2 . and Suspect 2.3.4 are very similar to those of SPheno 2.2.3. 
For illustration the results obtained with Suspect 2.3.4 are listed in Table |21 Here the 
scale dependence of mP-.f'^ is slightly larger, because the self-energy corrections are applied 

X2,3 

in the approximation given in Sect. 4.2 in [12]. 

The scale dependence gives an estimate of the size of the missing higher-order loga- 
rithmic corrections. This is not necessarily the full theoretical uncertainty. One point of 
caution is, for instance, the fact the SUSY threshold corrections to the gauge and Yukawa 
couplings are computed collectively at Mz- This is a valid choice. However, the matching 
between SM and MSSM couplings could also be done at, for instance, mLSP or Msusy- 
The thus induced uncertainty is not taken into account by the scale dependence as it 
results in a shift in the boundary conditions. 
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Case 


Ml 


M2 












X.-5 


(0) 


1 


(1) 


1 

X?. 


Xa 


A 


99.5 


192.4 


351.2 


96.2 


176.4 


358.0 


377.0 


95.2 


11 


iO.o 


357.0 


377.5 


B 


99.5 


192.4 


351.2 


96.2 


176.4 


358.0 


377.0 


97.9 


U 


^2.0 


357.6 


377.9 


C 


103.2 


192.9 


345.1 


99.9 


176.4 


351.5 


371.4 


101.5 


U 


^1.7 


351.6 


372.6 


D 


101.7 


192.1 


350.9 


98.4 


176.2 


357.9 


376.9 


97.3 


U 


^0.2 


356.7 


377.2 



Table 3: Isajet results (in GeV) for the neutralino sector at SPSla. Case A is the original 
Isaj et 7 . 72; Case B is Isaj et 7 . 72 with the improvement that the one- loop self-energies 
are each computed at their relevant scale as explained in the text; and Case C employs 
step-beta functions for all SUSY parameters. In case D, the SUSY parameters are all 
frozen out at Msusy = 456 GeV, and the one-loop corrections are applied at this scale. 

3.3 Step-beta functions in Isajet 

In the implementation employed in Isajet 7.72, each SSB parameter (aside from Higgs 
potential parameters) is extracted from the RG running at the scale equal to the SSB 
mass value. In particular, the parameters in the neutralino sector are Mi (Mi), M2(M2), 
/i(MsusY) and tan /9(Msusy) '^^ In this way logarithmic threshold corrections are included, 
and diagonalizing the tree-level mass matrix of Eq. (^J) gives a leading-log approximation 
of the neutralino pole masses. Finite corrections, however, can be of the same order as 
the logarithmic ones, so they have to be taken into account according to Eq. 

At this point some comments are in order on the actual implementation of the step- 
beta function approach in Isajet. Up to and including Isajet version 7.72: 

• While the SSB parameters are extracted from the RG running at their respective 
mass scale, the SSB parameters are not formally 'integrated out', so that the soft 
term RGEs remain those of the MSSM all the way from Mgut to Mz (unlike the 
case of the gauge and Yukawa couplings, where the beta functions change each time 
a threshold is passed). Thus, the Isajet algorithm is actually a mixed scheme. 

• For the finite shifts, the full expressions of [12] for the one-loop self-energies are 
used. These involve Aq, and Bx functions, which depend on the renormalization 
scale Q and are evaluated at Q = Msusy in Isajet 7 . 69-7 . 72. This leads to a 
double counting of logs between Msusy and the actual mass scale of the sparticle. 

Numerically Isajet 7.72 agrees quite well with the other public codes. The exception 
is the LSP mass, which, as already mentioned in the introduction, typically turns out 
to be a few per cent smaller than the results from SoftSusy2.0, Suspect 2.3.4 and 
SPheno 2.2.3. This can be understood as an effect of the above mentioned double count- 
ing, which induces an error oc ln(m^/M|usY)) where m is the mass parameter of the 
sparticle considered. Obviously, the effect is largest for the LSP. 

As a concrete example, the results of Isajet 7 . 72 for the parameters of SPSla, Eq. (jSj), 
are given as Case A in Table El Here m.'o* corresponds to the 'leading-log' approximation 

Xi 

while Wb\l is the final result including the 'finite' corrections. We note a difference of about 

Xi 

''inside Isajet, AfsuSY is called HIGFRZ and given by 
HIGFRZ=SQRT(MAX(AMZ**2 , AMTLSS*AMTRSS*SIGN(1 . , AMTLSS*AMTRSS) ) ) . 
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6 GeV in yu(MsusY) with respect to SPheno 2.2.3, which is however not of immediate 
concern for this analysis (it mainly stems from the different loop order of the effective 
Higgs potential in the two programs). More important for us is the 2 GeV difference in 
the LSP pole mass due to the choice of scale for the self-energy corrections. 

As a first step of improvement, we modify Isajet7.72 so that the one-loop self- 
energies for each neutralino (and also for each other sparticle) are computed with the 
renormalization scale set to its own mass scale. This means always setting the variable 
XLAM in the routine SSMILP to the mass of the sparticle being renormalized. In this 
way, the double counting is much reduced. The result is shown as Case B in Table El As 
expected, while the effect on the heavy neutralino masses is small, there is a shift upwards 
of 2 GeV in the LSP mass. 

Next we invoke the complete step-beta functions of Ref. [15] into the RGEs of all 
the SUSY parameters. The result of this is shown as Case C in Table El As can be 
seen, the agreement with the other codes is less good in this case. One reason is that 
when integrating out one parameter, this leads in principle also to finite shifts for the 
parameters that remain in the RGEs. This is, however, not taken into account in the 
'naive' step-beta function approach, which employs continuous matching conditions. We 
will discuss this and other issues in more detail in Sect. El 

As a consistency check, we also freeze out all SUSY parameters at Q = Msusy in 
Isajet, and compute the pole masses at this scale analogous to what is done in the 
other codes. This is shown as Case D in Table El and indeed gives the expected level of 
agreement with the other codes. 



4 Scalar top sector 

We next discuss the scalar top (stop) sector. At the tree level, the stop mass matrix is 
given by 

' \atht ml) y {^Atv<,-iiv^)ht M?^ + + J ' ^""^ 

where -Dj^ /j denote the D-term contributions, V\^2 are the Higgs VEVs and ht is the top 
Yukawa coupling. Using ht = mt/v2, the off-diagonal element can also be written as 

atht = {At- fi cot P) mt , (7) 

with mt the running top-quark mass. The difference to the neutralino sector is that 
this off-diagonal element can be large, actually as large or even larger than the diagonal 
elements. It can hence introduce a large mixing of tL and tu, and much enhance the 
splitting of the mass eigenstates ii^2- 

This is not a problem in the common scale approach, where all the parameters are 
taken at the same scale and therefore 

(_y^2)onshell _ y^2^g) ^ AA^?(Q) , (8) 

with AA^?(Q) the one-loop self-energy corrections given in [12]. Again, they have to be 
computed once for ti and once for ^2- This also results in a small (9(1%) difference in the 
stop mixing angle, depending on whether it is taken at = m? or = mf . 
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Q 


^Q3 




At 


DR 

m- 

ti 


DR 




pole 


pole 

t->. 


^pole 


100 


533.8 


454.6 


-524.4 


414.1 


607.0 


56.5° 


395.8 


577.0 


56.4° 


200 


516.3 


438.7 


-508.9 


398.7 


588.5 


56.4° 


399.6 


582.2 


56.4° 


^SUSY 


495.5 


419.9 


-490.4 


380.2 


566.2 


56.3° 


400.6 


586.0 


56.4° 


1000 


479.6 


405.5 


-476.3 


366.2 


549.6 


56.2° 


399.0 


585.2 


56.3° 



Table 4: SPheno 2.2.3 results (masses in GeV) for the scalar top sector at SPSla. 



Case 






At 


(0) 


(0) 

h 




(1) 
m- 


(1) 

to 


^(1) 


A 


493.7 


422.5 


-496.6 


381.2 


566.4 


55.7° 


401.8 


583.5 


56.8° 


B 


493.6 


422.4 


-496.4 


381.2 


566.4 


55.7° 


405.7 


584.9 


57.0° 


C 


488.5 


414.3 


-489.7 


373.7 


560.1 


56.0° 


394.1 


577.2 


57.2° 


D 


495.3 


420.6 


-498.5 


379.8 


566.6 


56.2° 


400.0 


583.9 


57.2° 



Table 5: Isajet results (masses in GeV) for the scalar top sector at SPSla. Case A is 
the result of Isajet 7.72; Case B is Isajet 7.72 with the improvement that the one- 
loop self-energies are each computed at their relevant scale; Case C employs step-beta 
functions for all SUSY parameters. In case D, the SUSY parameters are all frozen out at 
Q = ^nii^mi^ ~ 456 GeV, and the one-loop corrections are applied at this scale. 

In the step-beta function approach applied in Isajet 7.72, the lowest-order mass 
matrix given in terms of M?^(M?^), M?^(M?^), A(A), /x(Msusy), t^i,2(MsusY) and 
/ij(Msusy), and it becomes quite involved to define the corrections to the on-shell scheme 
in a consistent way. It is hence more convenient to integrate out all stop parameters at the 
average scale Qi = Msusy = and add the self-energy corrections as in Eq. (jH}. 

Table m shows the results of SPheno 2.2.3 for the stop sector at SPSla for different 
renormalization scales Q. As can be seen, the scale dependence of the pole masses at one 
loop is about 1.5%, that is 0{a1). 

TableElshows the results of Isajet for Cases A-D, analogous to Sect. 3.3. Case B has 
a somewhat higher ti mass than Cases A,D and SPheno 2.2.3. This can be explained by 
the fact that the ti and t2 masses are not respectively evaluated at = m?^ and = m?^ 
in Isajet 7.72, but both at = f^l^- Case C has somewhat lower top squark masses. 
The reason apparently is that the gluino has been integrated out at Q — 600 GeV, and 
its effect on the squark soft masses is to increase their values during running from Mqut 

5 Improved Isajet: effect on the complete spectrum 

In order to improve on the problems discussed above, the new version of Isajet, v7.73, 
adopts the following scheme: 

1. The running parameters of non-mixing sparticles, i.e. squarks and sleptons of the 
first and second generation, and the gluino, are extracted at their respective mass 
scale. The one-loop radiative corrections are implemented at this scale, which is 
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Mass 


Isaiet 7.72 


Isaiet 7.73 


SPheno 2.2.3 


^scale 


X? 

A. i 


95.19 


97.39 


97.11 


0.3 




180.5 


180.4 


180.7 


1.1 




356.7 


358.7 


364.9 


0.6 


xl 


377.2 


379.0 


382.2 


0.3 


xf 


180.4 


180.3 


180.3 


1.1 


xt 


376.2 


378.0 


383.3 


0.4 




203.2 


203.3 


202.4 


0.3 




144.0 


142.5 


144.1 


0.8 




187.1 


185.5 


186.2 


0.2 


fl 


134.8 


134.6 


134.4 


0.6 




206.7 


205.9 


206.4 


0.3 


Ur 


186.2 


183.4 


185.3 


0.2 


Ul 


559.5 


564.9 


565.1 


9.8 


Ur 


544.0 


548.6 


547.8 


8.9 


dL 


565.2 


570.9 


570.5 


9.7 


dR 


543.7 


548.2 


547.8 


8.9 


ii 


401.8 


395.2 


400.6 


5.0 


i2 


583.5 


584.4 


586.0 


9.0 


bi 


516.5 


518.6 


514.9 


7.9 


h 


539.7 


547.1 


547.5 


8.7 


9 


611.4 


605.9 


604.3 


1.3 



Table 6: Results of Isajet and SPheno for SPSla. Isajet7.73 is the new version, 
which includes the improvements explained in the text; S^'^^^^ is the scale dependence for 
Q = 0.1-1 TeV in SPheno 2.2.3. All values in GeV. 

also taken to be the renormalization scale. Thus, double counting of logarithmic 
corrections is not present. 

2. The parameters of mixing sparticles, i.e. neutralinos, charginos, stops, sbottoms, 
and staus, are all extracted at the scale Q = ^m^^m^^, which is also taken to be 
the renormalization scale. Diagonalization of the mass matrices are performed once 
for each sparticle mass, with self-energies evaluated at the corresponding sparticle's 
tree-level mass. 

3. The implementation of variable beta functions for the SUSY parameters is post- 
poned until a consistent treatment of logarithmic and finite corrections for multiple 
scales is available. 

In addition, in Isajet 7.73, the gluino mass radiative corrections depending on squark 
mixing have been added; these mixing corrections were absent in previous Isajet versions. 
The results of Isajet 7.73 are given in Table IHl and compared with those of Isajet 7.72 
and SPheno 2.2.3. Also shown is the scale dependence in SPheno 2.2.3. 

At this point it is important to note that the remaining differences in the spectra of 
Isajet 7 . 73 and SPheno 2.2.3 cannot be attributed exclusively to the different methods 
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GUT output 


Isajet 7 


72 


Isajet 7 


73 


SPheno 2.2.3 


A^QUT 


2.28 X 10^^ 


GeV 


2.28 X 10^^ 


GeV 


2.46 X 10^^ GeV 


gAMc) = g2{MG) 


0.715 




0.715 




0.721 




0.706 




0.706 




0.707 


htiMc) 


0.505 




0.516 




0.527 




0.049 




0.047 




0.051 




0.068 




0.068 




0.068 



Table 7: Comparison of DR GUT scale gauge and Yukawa couplings obtained from 
Isajet 7.72, Isajet 7 . 73 and SPheno 2 . 2 . 3 for the benchmark point SPSla. 

of implementing SUSY thresholds to the SUSY parameters. There are analogous differ- 
ences in the implementation of supersymmetric thresholds to the gauge and Yukawa cou- 
plings. In both programs, the experimental central values of the SM gauge couplings and 
third-generation fermion masses are used to determine the low-energy boundary condi- 
tions for the gauge and Yukawa couplings in the DR scheme. These are then extrapolated 
up to the GUT scale, defined as the scale where ai = 5ae/3(l — s^y) and a2 = Oie/s^/ 
unify. SPheno, as well as SoftSusy and Suspect, use a one-step implementation of the 
supersymmetric (log+finite) corrections to gauge and Yukawa couplings with the match- 
ing between SM and MSSM couplings done at the scale Mz- In Isajet, on the other 
hand, logarithmic thresholds are implemented in a one-by-one decoupling from the RG 
equations each time a threshold is passed, while finite corrections are implemented collec- 
tively at a common scale. The boundary conditions at the GUT scale are therefore not 
going to be identical. This is exemplified in Tabled where we show the GUT-scale values 
of the gauge and Yukawa couplings from Isajet and SPheno at SPSla. We see that there 
is no perfect gauge coupling unification since ^ ai = a2. The unification degree is of 
the order of 1.5%. This is expected for two- loop RGEs and can be attributed to threshold 
effects due to particles with GUT scale masses. The differences between Isajet 7.73 
and SPheno 2.2.3 amount to 8% for the GUT scale (which enters logarithmically in the 
RGEs), 1% for the unified gauge coupling, 2% for the top and 8% for the bottom Yukawa 
couplings. We conclude that even for a 'well behaved' point such as SPSla the differences 
in the GUT scale output cannot be neglected; in studies of Yukawa unified models they 
can become crucial. 



6 Discussion of the multiscale approach 

Integrating out all SUSY particles at a common scale is a reliable procedure if their 
masses are all in roughly the same ballpark. This is for example the case for the SPSla 
benchmark point. However, in the case where the SUSY spectrum [including the Higgs 
bosons] is spread over a large range of masses, one should get more precise predictions 
by integrating out the SUSY particles at various scales. Consider for instance the gluino 
mass. In the common scale approach 

mf^'CS=M3(Q) + Am5(g), (9) 
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while if the gluino is frozen out at its own mass scale, 



^poie,M. ^ ^^(^^^ ^ ^mg{M,) . (10) 

At SPSla with Mq ~ 550 GeV and M3 ~ 600 GcV, we find m?°''''^^ = 604.3 GeV 
and m~°^'^''^^ = 604.1 GeV, i.e. excellent agreement between the two methods. If we 
move, however, M3 to M3(M3) = 1.8 TeV, we find m?°'^'*"^ = 1512 GcV, and m^~°^^'^^ = 
1531 GeV — i.e. a spread of 19 GeV. A multiscale treatment of thresholds therefore seems 
desirable when the spectrum is considerably split. 

In the following, we outline the implications of a consistent multiscale approach. To 
work out a complete prescription is, however, beyond the scope of this Letter. Technically, 
at each scale a new effective theory (EFT) has to be constructed, and one faces the 
following difficulties: 



1. Finite shifts 

The naive way to take out particles from the RGEs via step functions (as done e.g. in [15]) 
and to use continuous matching conditions for the remaining parameters holds only for 
the lowest-order RGEs. At higher orders, finite shifts have to be introduced, as has been 
known for a long time [17-21]. The simplest example is gauge coupling unification in 
SU{5) theories: at lowest order, i.e. using one-loop RGEs, the boundary conditions at 
Mgut are given by g^^i) = gsu{2) = 5'su(3) = 5'su(5)- At next-to-leading and higher orders, 
finite threshold corrections due to particles with masses of order Mqut have to be taken 
into account [17,20,21], spoiling the lowest-order equahty of the gauge couphngs. Another 
prominent example is the evolution of the strong coupling between the scale of the lightest 
quarks and mz' see e.g. [22]. 

In our case, we are working with two-loop (or higher loop) RGEs for the SUSY param- 
eters, as for the gauge and Yukawa couplings. In this case, when integrating out SUSY 
particles at various scales, two issues have to be taken into account: 

(i) The shifts of the DR parameters of the sparticle that is integrated out to its pole-mass 
parameters involve contributions of all particles, which are degrees of freedom of the cur- 
rent EFT. 

(ii) The field(s) that are integrated out also lead to finite shifts in the boundary conditions 
of the parameters which remain in the RGEs. Take as an example the case where the 
gluino is the heaviest MSSM particle. At Q = |Af3| the gluino can be integrated out and 
the pole gluino mass rrig is obtained by taking into account the shifts from all strongly 
interacting particles (at the one-loop level; at the two-loop level all particles contribute 
in principle). In addition, there will be finite shifts for as, the Yukawa couplings (or 
equivalently the masses) of the quarks, and the squark parameters. A similar effect has 
been known in QCD for a long time, where the decoupling of heavy quarks also leads to 
shifts for the boundary conditions of running masses of the lighter quarks [23]. 



2. New couplings 

By integrating out part of the spectrum it may happen that the symmetry of the EFT is 
'smaller' than the symmetry of the underlying theory. In this case, there arise additional 
parameters in the EFT. Moreover, there are in general higher-dimensional operators com- 
patible with the reduced symmetry; well known examples of this are the neutrino mass 
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operator in the see-saw mechanism [24] or the operators governing rare meson/baryon 
decays. 

In the case of supersymmetric theories, there is a well studied and important example, 
namely the Higgs sector. Within the MSSM, supersymmetry requires that three of the 
seven couplings of a general 2-Higgs doublet model be zero and that the remaining four be 
expressed by gauge couplings (at least at tree level). Integrating out the SUSY particles 
does not only lead to finite shifts for the four non-zero couplings, but also introduces 
non-zero values for the couplings which are zero due to supersymmetry [25]. In the very 
same sector, a further complication arises if the top quark is integrated out, as SU(2) is 
broken in such a case. This leads to various new couplings to the W and Z bosons, each 
governed by its own RGE [25]. 

A second important example in this context is the bottom Yukawa coupling. If one 
integrates out, e.g. the gluino, from the spectrum, one induces non-holomorphic Yukawa 
couplings between the quark fields and the Higgs fields, e.g. a bbHu coupling [26]. 

The complete set of superfield operators leading to dimension 5 and 6 operators in the 
Lagrangian has been worked out for the MSSM [27], assuming that it is the effective theory 
of a more fundamental one. If one regards the MSSM as fundamental and integrates out 
part of it, the operators obtained will form a subset of those given in [27]. The techniques 
to obtain the tree- level coefficients in front of these operators have been elaborated in [28]. 



3. Sparticle mixing 

Clearly, once from a set of mixing particles the heavier ones are integrated out, the 
question arises of how to obtain the mixing effects that are visible in the case of a one- 
step decoupling. The answer is that the above mentioned higher-dimensional operators 
introduce these effects. Let us sketch this for the stops, for the case when m? ^ m? . 

Integrating out ti yields, for example, effective operators of the form 

C,i^tPRW,HltR , C2^^iPRWsH',iR (11) 

tL IL 

cJ-^H^uHTtRtl , cJ-^H',HTiRi*^ , C5^^i/M*t«tl + h.c.. (12) 

tL tL tL 

After electroweak symmetry breaking, the operators of the first line contribute to the 
titXi vertices, whereas those of the second line give contributions to the mass of the 
lighter stop. The coefficients Cj contain the information of the evolution of the operators 
from the scale m^^ down to m^^ . As a check that these operators really mimic the effect of 
mixing, let us calculate the stop mixing angle for this case, taking into account that the 
left stop mass is much larger than the m^^ and rrit (and thus the off-diagonal element): 

n_ -rritjAt - ficotp) -rritjAt - /icot/?) 

COS — y — — 2 • \ -L'J ) 

'm?^ — m|)2 + mf{At — fi cot f3y "^t^ 

At the scale m^^, Cj = 1 Vi, and thus if one naively replaces the Higgs fields by Vu,d/ 
and tfj by ti one obtains exactly the coupling of the lighter stop to the neutral wino in 
Eq. (fTTjl and the contribution of the left stop to the lighter stop mass in Eq. (fT^ . We note 
that the case of Cj = 1 is also obtained if one applies the see-saw formula to the stop sector. 
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7 Conclusions 



We have discussed the two general methods of treating threshold effects in the compu- 
tation of sparticle pole masses from high scale input: (a) renormalization at a common 
SUSY scale and (b) freezing out each sparticle at its own mass scale. We have focused 
in particular on the concrete implementations of these methods in SPheno 2.2.3 and 
Isajet 7 . 72 and compared the results of these programs for the SPSla benchmark point. 
Several shortcomings were pointed out, together with suggestions for improvements. The 
new version Isajet 7.73 incorporates these improvements, leading to better agreement 
between the various codes, especially for the LSP mass but also for the squark and gluino 
masses. 

Integrating out all SUSY particles at a common scale, as done in Sof tSusy, SPheno 
and Suspect, is a reliable procedure if their masses are all in roughly the same ballpark. 
The results of the alternative method followed in Isajet are in good agreement with this 
procedure. (Here, note that the method applied in Isajet is not a complete multiscale 
but a hybrid approach.) In order to match the ultra-high accuracy expected at the ILC, a 
consistent multiscale treatment of SUSY threshold effects seems necessary. The method 
is in principle well known from QCD. At each scale where a threshold is passed, the 
relevant field(s) have to be taken out of the RGEs and a new effective theory has to be 
constructed below the threshold. The shifts of the DR to the pole-mass parameters of the 
sparticle that is integrated out involve contributions of all the particles that are degrees 
of freedom of the current EFT. Moreover, the field(s) that are integrated out lead to 
finite shifts in the boundary conditions of the parameters which remain in the EFT, so we 
face non-trivial matching conditions. Additional complications arise when the symmetry 
of the EFT is 'smaller' than the symmetry of the underlying theory. We have outlined 
this multiscale approach in Sect. 6 and discussed its technical implications. A complete 
prescription for the matching over multiple scales is, however, beyond the scope of this 
Letter and will be explored in the future. 
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